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Aims. A solution of the radiative-transfer problem in 3D with arbitrary velocity fields in the Eulerian frame is presented. The method 
is implemented in our 3D radiative transfer framework and used in the PHOENIX/3D code. It is tested by comparison to our well- 
tested ID co-moving frame radiative transfer code, where the treatment of a monotonic velocity field is implemented in the Lagrangian 
frame. The Eulerian formulation does not need much additional memory and is useable on state-of-the-art computers, even large-scale 
applications with 1000's of wavelength points are feasible. 

Methods. In the Eulerian formulation of the problem, the photon is seen by the atom at a Doppler- shifted wavelength depending on 
its propagation direction, which leads to a Doppler- shifted absorption and emission. This leads to a different source function and a 
different A* operator in the radiative transfer equations compared to the static case. 

Results. The results of the Eulerian 3D spherical calculations are compared to our well-tested ID Lagrangian spherical calculations, 
the agreement is, up to v max - 1 • 10 3 km s _1 very good. Test calculation in other geometries are also shown. 
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1. Introduction 

A solution of the ID radiative transfer problem in arbitrary 
vel ocity fields in the L agrangia n frame has been developed 
by |Baron & Hauschildt| (2004) . |Knop et al.] ( |2009| ) optimized 
their method by reducing the memory footprint of the algorithm 
(using domain decomposition), and also introduced a new 
method which speeds up the formal solution by developing 
an iterative Gauss-Seidel (GS) type solver where the solution 
becomes quasi-analytic when the source function is interpolated 
linearly. However, in a 3D setup the limiting factor is the 
memory footprint which already stretches the limits of modern 
supercomputers. The memory requirements of a 3D calculation 
in the Lagrangian frame are very high because one has to store 
additional wavelength information in every volume element 
(hereafter voxel) for every solid angle, whereas the 3D Eulerian 
frame calculation needs more computing time in multi-level 
applications due to the explicit computing of the opacity for 
every solid angle point. 

In Section [2] we describe the Eulerian formulation of the 
problem, in Section[3]the comparison to our well-tested ID code 
is presented and application examples in other geometrical se- 
tups are shown. In Section [4] we describe expected and discov- 
ered limitations with the formalism in the Eulerian frame. 



straight lines with given direction (#, 0) in the Eulerian frame, 
is simply given by 
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With this definition, the formal solution for a given characteristic 
(#, 0) can be written as 
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where i labels the points along the characteristic, S is the source 
function and Ar t is the optical depth, computed, e.g., by using 
piecewise linear interpolation of the opacity x along the charac- 
teristic: 



Atz-i = (Xi-iW +XiW)\st-i ~ st\/2. 



(4) 



Hauschildt & Baron (2006) give a more detailed explanation of 



the general method. 

Due to the movement of the atom in the Eulerian frame, the 
atom 'sees' the photon (on the characteristic) at a wavelength 
shifted according to 
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2. Method 

The 3D radiative transfer framework uses the full-characteristics 



method to solve the radiative transfer equations (Hauschildt & 



where e C har is the unit vector in the direction of the characteristic, 
c the speed of light and v the velocity of the atom. This leads to 
a different opacity seen by the characteristic depending on its 
direction: 



Baron| |2006| ). The intensity along a characteristic, which are XiW =Xi(A,0,(f>) 



(6) 
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In the case of line transfer, the profile of the line becomes 
anisotropic in the Eulerian frame: 



O(/t) = OOl,0,( 



(7) 



The solution for the line transfer then proceeds with the Eulerian 
profile function 0(/l, 6, 0) and following |Baron & Hauschildt 
420071 ), we obtain 



/ = III 0(^, fi, (p)J A dAdjid(p 
and 
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where A is the wavelength in the observers frame, \i - cos 6 and 
the solid angle under which the voxel is hit by the characteris- 
tic. 

Relativistic corrections are neglected in the Eulerian frame, as 
we are working on a implementation of the co-moving frame 
method into our 3D code ( [Chen et alT] |2007| E. Baron et al, in trum from the ID code, 
preparation). 



3. Application Examples 

A simple two-level-atom approach was used to test the code, the 
physical atmosph ere setup presented in Sectio n [3.1|3.2 is similar 
to the one used in Hauschildt & Baron (2009), remarks about the 
numerical resolution of our spherical 3D code can also be found 
there. 



Fig. 1. Comparison of the Eulerian 3D and the Lagrangian ID 
code with a linear increasing velocity field with v max = 1 • 
10 3 km s" 1 and scattering 6ii ne = 10" 2 . Solid: Spectra from all out- 
ermost voxels in the 3D spherical coordinate grid, dashed: spec- 



3.1. 3D Spherical coordinates 

Our well-tested Lagrangian ID spherical code was used to com- 
pare its results to the new Eulerian formalisms in the 3D RT code 
to verify the 3D Eulerian Code. 

For this, many test calculations were made with linearly in- 
creasing velocity fields up to v max = 1 • 10 3 kms -1 . Figure [l] 
shows the corresponding spectra from a calculation with v max = 
1 • 10 3 kms -1 , the agreement to the ID Lagrangian code is 
very good, the variation of the 3D lines is due to numerical 
resolution. The 3D model was computed using (n^ng,^) = 
(197x99x197) = 3.842.091 voxels and a solid angle resolution 
of (0 C , C ) = (64, 64) (computing time: ca. 18h on 2048 CPUs)[] 

3.2. 3D Cartesian coordinates with and without Periodic 
Boundary Conditions 

In Figure [2] a test calculation of the 3D Cartesian code with pe- 
riodic boundary conditions in a (n x ,n y ,n z ) = (157x157x157) = 
3.869.893 voxel grid with (0 C ,</> C ) = (64,64) is shown (com- 
puting time: ca. 6h on 2048 CPUs) l . The difference between 
the calculation with and without velocity field is clearly visible. 
Similar results can be obtained in the Cartesian mode without 
periodic boundary conditions, the spectra are not shown here. 



3.3. 3D Hydro-structure 

To test the new method with a 3D hydro-dynamical structure 
with an inherent arbitrary velocity field, we obtained a computed 
snapshot of convection in the solar atmosphere from Ludwig 
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Fig. 2. Spectra from all outermost voxels in the 3D Cartesian 
geometry with periodic boundary conditions, a linearly increas- 
ing velocity field with v max = 1 • 10 3 km s" 1 and line scattering 
£iine = 10~ 2 - Solid: Without the treatment of the velocity field, 
dashed: with the treatment of the velocity field. This plot shows 
the influence of the velocity field on the line. 



( |Caffau et alT] [2007} |Wedemeyer"et~aL| [2004] ), which was used 
as input for the code. The spectra from a few outermost vox- 
els in a calculation ((n x ,n y ,n z ) = (140x150x140) = 2.940.000, 
(0 c ,0c) = (64,64), computing time: ca. llh on 1024 CPUs l ) 
with and without velocity field is shown in Figure [3] The plot 
show that the use of velocity fields in such calculations is manda- 
tory. 

4. Limitations of the Eulerian formalisms 

4.1. Solid Angle points 

In static line transfer problems it is necessary that the profile of 
the line is covered by the discretized wavelength grid used in the 
calculation, as the profile does only depend on the wavelength 
O = O(zl). In Eulerian moving atmospheres the profile depends 
on equation [7] and therefore also on the solid angle discretiza- 
tion. 
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Fig. 3. Spectra from selected voxels on top of the 3D hydro-dynamical Cartesian atmosphere to show the influence of the inherent 
velocity field on the line, the calculations include line scattering (e^ = 10" 2 ). Solid: without velocity field, dashed: with the inherent 
velocity field. 




4.2. Relativistic Velocities 

The formulation of the Eulerian method in the observers frame 
is inherently non-relativistic, which leads to differences between 
the ID Lagrangian code and the new method presented here 
when the velocity field is greater then about 5 • 10 3 km s" 1 . When 
the velocity is getting close to the speed of light, the Lorentz 
boost drives the continuum of the radiation field higher, what is 
clearly visible in Figure [5] As v gets closer to c this effect in- 
creases, the error in the Eulerian solutions increases rapidly. 




Fig. 4. Plot to illustrate the problem of a poor solid angle dis- 
cretization. The plus signs show the profile in the line center 
of a static atmosphere. The asterisks and the diamonds show a 
poor/medium solid angle discretization, the triangles good solid 
angle discretization in the Eulerian moving atmosphere where the 
anisotropic Eulerian profile has good coverage. See text for a more 
detailed explanation. 
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In Figure [4] the profile in the (observer's frame) line center is 
plotted: The plus signs show the profile in the static case, as it is 
not solid angle depend it is everywhere 1 . The asterisks show the 
anisotropic profile of a poor quality, the diamonds of a medium 
quality solid angle discretization in the moving Eulerian atmo- 
sphere: The profile in the line center is not hit at all, this causes 
a wiggly or even no spectral line. The triangles show a good 
quality solid angle discretization, the anisotropic profile in the 
Eulerian moving atmosphere hits the line center and the profile 
is well covered. 

Various test show that a solid angle resolution of (0 C ,0 C ) = 
(64, 64) is sufficient for velocity fields up to v max = 10 3 km s _1 . 



Fig. 5. Example of a calculation with a linearly increasing velocity 
field with v max = 1 • 10 4 km s" 1 . Solid: Spectra from all outermost 
voxels in the 3D spherical grid, dashed: Spectra from the ID code. 
The Lorentz boost drives the continuum up as v gets closer to c. 



Various tests show that the agreement between the well- 
tested Lagrangian ID and the Eulerian 3D code is excellent with 
velocity fields up to 1 • 10 3 km s~ , this is enough to do 3D radia- 
tive transfer in convection or global circulation models. We have 
extended the method described in Knop et aL] ([2()09) into our 3D 
code, which then allows exact, full relativistic radiative transfer 



in 3D ( Chen et al. 2007 E. Baron et al., in preparation) 
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5. Conclusion 

With our new Eulerian method it is now possible to do 3D radia- 
tive transfer in non-relativistic, arbitrary velocity fields in spher- 
ical, Cartesian, and cylindrical (while not described in detail, the 
method is also implemented in this part of the code) coordinates. 
The memory footprint and the computing time of the new algo- 
rithm in the two-level-atom setup presented here is negligible 
compared to the general requirements of the 3D code. In multi- 
level applications the time needed to calculate the opacity for 
every solid angle must be considered. The velocity field limita- 
tions must also be kept in mind when using the method. 
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